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^ ■ Abstract 

A theory of gravitational quantum states of ultracold neutrons in waveguides with absorb- 



ed 



■ ing/scattering walls is presented. The theory covers recent experiments in which the ultracold 
neutrons were beamed between a mirror and a rough scatterer /absorber. The analysis is based 

J> ' on a recently developed theory of quantum transport along random rough walls which is modified 

m 

' in order to include leaky (absorbing) interfaces and, more importantly, the low-amplitude high- 

CN ; 

■ aperture roughness. The calculations are focused on a regime when the direct transitions into the 

O 



continuous spectrum above the absorption threshold dominate the depletion of neutrons from the 
gravitational states and are more efficient than the processes involving the intermediate states. 
The theoretical results for the neutron count are sensitive to the correlation radius (lateral size) 
of surface inhomogeneities and to the ratio of the particle energy to the absorption threshold in a 
weak roughness limit. The main impediment for observation of the higher gravitational states is 
the "overhang" of the particle wave functions which can be overcome only by use scatterers with 



> 

■ strong roughness. In general, the strong roughness with high amplitude is preferable if one wants 

just to detect the individual gravitational states, while the strong roughness experiments with 
small amplitude and high aperture are preferable for the quantitative analysis of the data. We also 
discuss the ways to further improve the accuracy of calculations and to optimize the experimental 
regime. 
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I. INTRODUCTION 



One of the recent discoveries in neutron physics was the experimental observation of quan- 
tization of motion of ultracold neutrons in gravitational field Q . By itself, the quantization 
of motion of particles by a linear potential is a well-known textbook problem in quantum 
mechanics and has already been observed for other ultracold particles (for example, for spin- 
polarized atomic hydrogen in magnetic field with a linear gradient However, the Earth 
gravitational field is so weak and the energies of the corresponding discrete quantum states 
for neutrons are so low (of the scale of 1 peV!) that the observation of such states is indeed 
a major experimental challenge. What is more, the measurements in this energy range open 
the door, for instance, to experimental studies of the fundamental short-range forces for 
neutrons near solid surfaces (see 0, |^ and references therein) . This unique opportunity is 
inherent to the high sensitivity of the neutron wave functions in the gravitationally bound 
quantum states of neutrons to the existence of any new type of short-range forces between 
a neutron and a surface. 

The schematics of a typical experiment is the following. A collimated beam of ultracold 
neutrons is sent between two, usually horizontal, solid plates at a distance several fim from 
each other. The material of the plates is such that the neutrons hitting the plates with the 
vertical velocity above, approximately, 4 m/s, get absorbed by the plates, while the particles 
with a lower vertical component of velocity get refiected. One of the plates, usually the upper 
one, has a rough surface while the other one is an almost ideal mirror. The quality of the 
mirror is surprisingly good: the single-refiection losses are in the range 10^^ -r- 10"^. As a 
result of scattering by a rough wall, the direction of the neutron velocity changes and the 
neutrons can acquire a sufficiently large vertical component of velocity to get absorbed by 
the walls. Only the particles, which have extremely small vertical velocity (low gravitational 
energy) do not reach the rough scatterer (the upper plate) and, therefore, are not scattered 
and absorbed. These particles get through the system and get counted by the exit counter. 
[In practice, the bottom mirror is often made longer than the scatterer/absorber and the 
entrance to the neutron counter is collimated with a long horizontal narrow slit. With such 
configuration, the absorption is not the only source of the neutron loss. The purpose of 
this slit is to avoid counting of neutrons which exit at a large angle to the mirror surface 
and could otherwise get counted, albeit with a small probability, either directly or after 
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scattering by spectrometer walls]. 

A quantum description of experiment differs from the classical one above in several im- 
portant aspects. The motion of neutrons in the gravitational field between the plates is 
quantized. The lower the energy, the smaller the (vertical) size of the state. The particles in 
the lowest gravitational states do not reach the scatterer (or, more precisely, reach it with an 
exponentially low probability), do not scatter, and, therefore, do not get absorbed. With de- 
creasing distance between the plates, the lower and lower states will reach the scatterer and 
get absorbed. As a result, the neutron count in the case of spatial quantization of neutrons, 
should, in ideal circumstances, be a step-wise function of the clearance between the plates. 
The positions of the steps should then provide the calibration for ultralow energy measure- 
ments. In experiment, unfortunately, such a clear-cut step-wise picture does not emerge: 
though the curves, without doubts, demonstrate the quantization of the gravitational states, 
the higher gravitational states have, so far, eluded identification. 

One of the main purposes of this paper is to develop a general theory that would cover the 
experiment, explain the difficulties associated with identification of the higher states, and 
suggest the ways to optimize the experimental regime. Many aspects of the problem have 
been discussed very recently in Ref. ^] which analyzed potential mechanisms of neutron 
losses in experiment Q|, including the surface roughness as the principal mechanism of non- 
specularity of refiection which is ultimately responsible for the neutron depletion. This 
paper demonstrated that the losses of neutrons from the lowest gravitational states are 
determined by the neutron tunneling through the gravitational barrier between the classical 
turning point and the position of the scatterer/absorber. 

However, the theory ^] did not provide an explicit quantitative link between the rate 
of neutron loss and parameters of the roughness of the neutron waveguide. Since the non- 
specularity of refiection by the rough surfaces is an obvious principal source of the neutron 
depletion in experiment jj , any accurate quantitative description of the phenomenon should 
express the neutron count directly via the parameters of the random surface roughness. Our 
results below not only provide such an explicit link, but also suggest what kind of roughness 
is the most advantageous for experiment. 

In order to extract quantitative information from the experiment and, eventually, to study 
the short-range forces, one needs an accurate theoretical description of quantum diffusion of 
neutrons along rough walls. Recently we developed a theory of quantum transport in quasi- 
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2D systems with random rough boundaries (see Refs. jaLD] and references therein). In this 
paper we apply this formahsm to quantized neutrons in the gravitational field. The practical 
purpose is twofold: to develop an accurate quantitative theory of neutron transport in the 
presence of random surfaces, which can be used for a quantitative analysis of experimental 
data, and to determine the parameters of surface roughness for the scatterer surface that 
would help to optimize future experiments. As for the purely theoretical benefits, we are 
able to extend our general theory to include strong, high-aperture roughness and to develop 
a theoretical description for systems with absorbing, "leaky" walls. 

Roughness of the walls affects the quantized, in this case, gravitational, states of particles 
inside the quantum wells by shifting and broadening the energy levels. The stronger the 
roughness, the larger are these shifting and broadening. In this paper we, effectively, cal- 
culate the line broadening which is determined by the roughness- driven transition between 
the states. Our perturbative approach works as long as the resulting uncertainty in the level 
positions is smaller than the distance between the levels. The same condition is necessary 
for the quantitative interpretation of experiment: though strong roughness can make the 
effect of the gravitational quantization of states more pronounced, the uncertainty in the 
level positions can render the results unusable for the calibration of the short-range forces. 

The strength of roughness is characterized by the ratio of its amplitude i to the size 
of the well L and by the aperture of roughness i/R, which is determined by the ratio 
of the amplitude of roughness i to its typical lateral scale R (the correlation radius of 
surface inhomogeneities) . Usually, both the line shifts and broadening are small as long as 
i/ L, i/ R 1 (see, for example, Ref. ^ and references therein). 

In this particular case there could be two exceptions. First, the highest energy levels, 
which are close to the absorption threshold, could become relatively broad even for a mod- 
erate roughness. This is not very important for the results below since in this paper we deal 
mostly with direct roughness-driven transitions from the gravitational states into the con- 
tinuous spectrum above the threshold and can disregard the transitions via the intermediate 
states. Second, and more important, the effect of surface roughness on the lowest gravita- 
tional states is anomalously small because for such states the particle wave functions on the 
rough walls are exponentially small in the distance between the classical turning point and 
the wall. Therefore, as long as we can disregard the absorption processes which involve the 
intermediate, higher levels, the linewidths of the lowest gravitational states remain small 
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and the calculation of the absorption rate can even be extended to stronger, high-aperture 
roughness. We plan to address the neutron absorption that involves the transitions via the 
intermediate states later. 

When interpreting the experimental results, the main uncertainty is introduced by the 
lack of accurate data on parameters of surface roughness, especially on the type of corre- 
lations and the correlation length of surface inhomogeneities. Simple measurements of the 
average amplitude and the lateral size of roughness are insufficient when one needs to know 
the shape of the correlation function of surface inhomogeneities which can be extracted, 
almost exclusively, from diffraction measurements j^. Another major sort of uncertainty is 
a relatively poor information on the energy/velocity distribution of neutrons in the beam. 
While the dispersion in horizontal component of velocity can be measured very accurately, 
the precise measurements of the vertical dispersion are much more problematic. Therefore, 
one of the goals of this paper is to identify the range of parameters in which this lack of 
precise input data affects the results as little as possible. Surprisingly, this is possible. 

The paper has the following structure. In the next Section we describe the neutron 
wave functions and eigenstates, introduce proper dimensionless variables, and present a 
set of transport equations for quantized neutrons near random rough wall(s). In Section 
III we present the roughness-driven scattering probabilities in the case of weak roughness. 
In Section IV we explore the ways to generalize our theory to stronger roughness with 
some success in the case of high-aperture low-amplitude roughness. The final results for 
absorption times for neutrons on the leaky walls and for the neutron count for direct and 
inverse geometries are given in Section V. Section VI contains conclusions and discussion of 
how to improve our theory and optimize the experiment. 

II. MAIN EQUATIONS 

A. Potential well, notations, and dimensionless variables 

In the bulk of the paper we will disregard the ultrashort-range potential for neutrons 
near solid surfaces and consider surfaces of solid plates as potential barriers for neutrons 
of the height Uc (the absorption barrier). We will discuss the short-range forces near the 
walls towards the end of the paper. If one forgets for a moment about the roughness of 
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FIG. 1: Schematic drawing of the potential well in vertical direction z. Dashed lines - the "real" 
potential; solid lines - the potential used in the paper. Since mgH/Uc ~ 10^^ and we are particu- 
larly interested in the lowest energy levels, the corrections are negligible [with an accurately scaled 
potential, the difference between the dashed and solid levels cannot be seen]. 

the scatterer/absorber (the upper plate), the potential energy for neutrons in the vertical 
direction, U {z) , is ?7 < 0) = Uc + mgz (mirror at the bottom), U (z) = mgz at < z < H 
(the gravitational field between the plates; the scatterer/absorber is at the distance H from 
the bottom mirror), and U (z) — mgz + Uc = Uh + fng {z — H) (the absorption barrier) 

z > H. Schematically, this potential is presented in Figure 1 by the dashed lines. An 
alternative geometry is discussed in Section VE. 

Numerically, Uc ~ 1. 34 x 10~^^ J (the kinetic barrier for absorption is 4 m/s). The 
characteristic gravitational energies inside the well, on the other hand, are by five orders 
of magnitude lower, mgH ~ 1.6 x 10~^^ J (for ~ 10 /xm). This means that the plot 
in Figure 1 is out of scale. Since we are interested almost exclusively in the gravitational 
states near the bottom of the well, we can, without any noticeable loss of accuracy, ignore 
the gravitational field near the top of the well and replace the exact potential (the dashed 
lines) by the one given in Figure 1 by solid lines. 
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The energy spectrum of neutrons in this potential well consists of discrete energy states e„ 
below Uc and continuum above. [Note, that the short-lived discrete states near the threshold 
could be also treated as continuum. In the well walls are not ideal, but rough, all the energy 
levels, especially the higher ones, broaden. With an increase in roughness amplitude, lower 
and lower levels will acquire a noticeable width. This is, in some sense, equivalent to the 
lowering of the absorption threshold with increasing roughness]. We will call the discrete 
states close to the bottom of the well the "gravitational states" and higher states in the 
main part of the well - the " square well states" . Though this terminological shortcut is not 
well-defined, it is convenient for discussing the physics and the results because of a very 
large ratio Uc/mgH ~ 10^; in computations, unless stated otherwise, all the states will, of 
course, be treated accurately. 

It is convenient to measure the distances z and the energies e„ in units of Iq and cq, 
An = Cn/eo and s = z/Iq, where /q = ^^^^ {2m?g) ~ 5.871 /im and Cq = frtglQ ~ 0.602 
peV~ 9.6366 x 10^'^^ J are the size and the gravitational energy of a neutron in the lowest 
quantum state in the infinite gravitational trap. In these units, the overall kinetic energy 
of particles in the beam E and the barrier height Uc are 1.4 x 10^ <e = -E'/eo<8.7x 10^, 
Uc = Uc/cq ~ 1.4 X 10^. A very important parameter is the ratio x = Uc/e. In experiment, 

0.15<x = uc/e<l. (1) 

The value of this parameter shows to what extent the kinetic energy of a neutron exceeds 
the barrier and determines how easy is it for a neutron to get absorbed by a plate when 
the direction of velocity is rotated as a result of scattering by the rough scatterer/absorber. 
As we will see, the results are very sensitive to the exact value of x even within the limited 
range (Q) because both values of Uc and e are very large. 

It is also convenient to introduce the dimensionless velocities (momenta) in the beam 
direction along the well (x direction) Pj, Vj = PjVo, Vq = y/2gl^ ~ 1.073 x 10^^ m/s, 
Pj = — Xj = Pjlo ~ (l — l-^i/^) = a/c — A j (for lower levels, Xj <^ e and ~ 
^/e (l — |Aj/e)). The range of kinetic energies in experiment is 373 < y/e < 932. 

The roughness of the surface of the scatterer/absorber (the upper plate) is characterized 
by the correlation function of surface inhomogeneities which is determined by the average 
lateral size (correlation radius) R and the average amplitude (height) £ of the inhomogeneities 
(for precise definitions see Section III). In dimensionless variables, r = R/Iq, r] = E/Iq. It is 
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reasonable to expect that 77 < 0.1 [so far, in experiments [l| 0.03 ^rj < 0.1; roughness with 
higher amphtude than 0.1 would be too strong for a quantitative theory]. We will discuss 
how to introduce the surface roughness into the equations later on. 

B. Wave functions and energies of the bound states 

At E < Uc (or A„ < Uc = Uc/cq) the wave function ipn (s), s = z/Iq, inside the well, 
0<s<h = H/lo, is 

(s) = A„ ■ [Ai {s - A„) - S-^Bi {s - A„)] . (2) 

Outside the well, 

'(pn (s) = Bn ■ exp [-kn {s - k)] , s> h, (3) 



Ipn (s) = Cn ■ exp [knS] , S < 0, kn = a/Mc " A„. 

The values of parameters Sn and the energy eigenvalues A„ are determined by the continuity 
conditions at s = and s = h: 

Cn = An ■ [Ai (-A„) - S„Bi (-A„)] , Bn = An- [Ai {h - A„) - S„Bi {h - \n)\ , (4) 

Cnkn = An ■ [Ai' (-A„) - 5„Bi' (-A„)] , Bnkn = -An ■ [Ai' {h - A„) - 5„Bi' {h - Xn)] , 

which are equivalent to two equations 

1/kn = [Ai (-A„) - 5„Bi (-A„)] / [Ai' (-A„) - 5„Bi' (-A„)] (5) 

= _ [Ai {h - Xn) - SnBi {h - Xn)] / [Ai' {h - Xn) " 5„Bi' {h - A„)] . 

Computation of the neutron count would also require the normalization coefficients An, 

A-^ = ]-f [Ai {s - Xn) - 5„Bi {s - Xn)f ds (6) 
+ [[Ai {-Xn) - SnBi {-Xn)f + [Ai {h - A„) - ^„Bi {h - A„)]'] , 

and Bni Cn-, Eq. (j3]). Some information on discrete states can be also found in Ref. 

C. Approximate description of the bound states 



a. 



For what follows, we will not always need the exact eigenvalues and the wave functions. 
First, let us note that the threshold energy Uc is much higher (by about five orders of 
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magnitude!) than the gravitational energy mgH in a typical experiment, Uc ^ h. One can 
also separate the states into two groups: the lowest states for which the presence of the 
gravitational field near the bottom of the well is important (the "gravitational states") and 
the higher states for which the fact that the gravity determines the shape of the bottom of 
the well is mostly irrelevant (the "square well" states). 

For the gravitational states in an infinite well, Uc — > oo, the wave functions Q should 
become equal to zero on the walls s = 0; /i, and equations for the eigenvalues Eq. (jSJ should 
be replaced by the following equation for the energy spectrum A„ {h) [here and below we 
identify physical parameters, which are calculated for a well with infinite walls, by a bar 
over the corresponding symbols]: 

Sn = Ai (-A„) /Bi (-A„) , Ai (/i - A„) - :S„Bi {h - A„) = 0. (7) 

The normalizing coefficients An are now defined as 

= j-«n, On = 1^ [Ai (s - A„) - S^Bi (s - A„) ] ^ rfs | (8) 

= I (Ai' (-A„) - 5„Bi' (-A„))' - (Ai' [h - A„) - 5„Bi' (/^ - A„))'}"\ 

If for the lowest gravitational states the size of the state A„ is much smaller than the spacing 
between the walls, A„ <^ h, the equation for the energy spectrum is even simpler. 



Ai -A, 



0. 



(9) 



with the semiclassical eigenvalues 



Af = — 2n - 



2/3 



and the normalizing coefficients 



Ai = —ar 



In 



Ails - Xr 



ds 



iAi'-2 

'0 



-A. 



(10) 



Since Uc is very large but finite, Uc ~ 10^, the eigenvalues A differ from their values A (|7j) 
at Uc ^ cxD, A = A + 5A, S = S + 5S . The corresponding corrections can be calculated by 
expanding Eqs. © in 1/ a/uc — A: 

1 Bi(-A) 



\/uc- \ 
1 



5\ + 5S 
5\ + 5S 



Ai' (-A) - 5Bi' (-A) ' 

Bi {h - A) 
Ai' {h-\)- SBi' (/i - A) ■ 



(11) 
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Then the normahzing coefficients B,C [B = C = 0), become 

B = ^Ij{H) = -A [Ai' {h-X}- SBi' {h-X}] /Ju,-X, (12) 



C = ^ (0) = A [Ai' (-A) - SBi' (-A)] /^Uc - A. 

For the square well states well above the bottom, but still below Uc, one can neglect the 
gravitational field near the bottom and use the standard equation for the eigenvalues A„, 

tan(VXh)=-y^^±^, (13) 

V / k-Vxs ^ 



k = ^u, - A, S = VA/ (uc-X). (14) 
For these symmetric states the values of the wave functions on the walls are 

(0) = (H) = A'S\ (15) 



-2 _ ^0 r,2 I ^0 



2 (^^/Xh + S + y/XhSy - 25 cos (2VXh^ - (l - S^) sin (2y/Xh 



k 4:^/X 

In the case of very deep levels X/uc, Eqs. ffT^ — (fT^ can be simplified as 

A = Afl--^y S=J-, (16) 



€{^) = €{H) = ^-, (17) 

where A and A are the eigenvalues and the normalizing coefficient for the infinitely deep 
well, 

A„ = Tx^n'/h^ aI = 1/ {2hh) . (18) 

Approximately, the total number of discrete levels is 

h 

M 0i;^12O/i. (19) 

TT 

D. Time evolution of the neutron population 

Recently we developed a versatile formalism for description of 2D diffusion of quan- 
tized particles along random rough walls j^. The key is to include the information on the 
wall roughness not into the boundary conditions, but into the roughness-driven transition 
probabilities between the quantum states in the transport equation. The application of 
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this formalism to the neutron problem at hand requires several important modifications: 
we should include neutron absorption by the "leaky" walls and reconsider the time/space 
properties of the transport equation. 

The possibility of neutron absorption by the walls can be taken into account by the 
calculation of the probability of transitions into the continuous spectrum without the coun- 
terbalancing terms in the transport equation. The second anomaly of this problem is that we 
are calculating not a stationary diffusion flow but the time dynamics of the neutron count. 
As a result, we have to solve a time-dependent problem with an initial condition in the form 
of initial distribution of neutrons over the vertical velocities, i.e., over the gravitational 
quantum states. 

The main modification deals with the dimensionality of the problem. In we considered 
the scattering-driven 2D diffusion in two equivalent directions parallel to the walls while the 
third dimension (the finite motion perpendicular to the plates) gave rise to a set of discrete 
quantum states. In our neutron problem there is a preferred direction along the plates 
(the beam direction x); the neutron count is being done only in this direction. Therefore, 
we should exclude the sideway motion and reformulate the transport problem as a 2D 
problem with vertical motion (coordinate z) and the direction of the beam (x-direction 
with the corresponding momentum p). In other words, we have to modify the transport 

quasi-2D particles (particles 

with continuous spectrum in x, y directions and discrete quantized gravitational states in z 
direction) in order to eliminate the momentum in q direction y. 

According to the distribution functions rij {q,p) in each quantized state j obey the 
following Boltzmann equations: 



(20) 



where the 2D vector q in the plane of the mirror has y and x components q = {q,p). Since 
we disregard the change of momentum in ?/-direction, 

27r 

Wjj' (q, q) = S{q- q) Wjj, {p,p') , rij (g, p) = —6 (g) rij (p) , (21) 
and the integration over dq' in Eq. (l2(Jj) yields 



27r ^ /" {p,p') {p') - n, {p)] 6 (p'V2m - pV2m + e,, - e,) , 
J' 
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(22) 



We are dealing with particles with a fixed energy E, i.e., with particles in gravitational 
quantum states ej with distinct values of the lateral momentum pj = ^y2m [E — Cj), 

n,{p) = —5{p-p,)N,, (23) 

where Nj is the number of particles in state j per unit length of the beam. Finally, the 
equations for the particle balance (j22|) reduce to 



^ = E ^n' iPjiPr) [Nf/vr - N,/v,] . (24) 

To this equation we should add the terms responsible for the neutron absorption (transitions 
into the continuous spectrum): 

dNj/dt= Wjy{Nj,/vy - Nj/vj) - Nj/vj f Wx'jdX'. (25) 

Note, that the separation of the energy spectrum into discrete lines and continuum, which 
is exact for an ideal well in Figure 1, is unambiguous only for weak roughness. In general, 
roughness leads to shifting and broadening of all the lines. The transition probabilities 
and, therefore, the broadening are expected to be higher for the higher levels (see the next 
Section). With increasing amplitude of roughness, the higher levels can become broad 
enough to be considered continual, which, in turn, is equivalent to the lowering of the 
absorption threshold. In the case of the high-amplitude roughness this effect could even 
become dominant. However, in the case of high amplitude roughness the whole energetics 
of the systems changes so dramatically that it could be impossible to use the experimental 
results in a meaningful quantitative way. 

The calculation of the neutron count can be done either via the time dependence of the 
neutron number at fixed position during the time of fiight t = L/v^ or via the decrease of 
the number of neutrons while they move along the absorber of the length L = tv^- Within 
the former approach, if initially there was Nj (0) neutrons per unit area in each discrete 
state, the time evolution of the densities Nj is described by the particle balance equations 

(1221). 

The set of linear differential equations describes the exponential disappearance of 
neutron with a set of characteristic relaxation times r. These relaxation times r are given. 
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according to Eq. (j2SI) , by the following characteristic equation: 



= det 




(26) 

' ^3 Jy>Uh ^33' ^3' ^3 ju^z 



where it is convenient to introduce the dimensionless velocities j3j, Vj = PjVq, Vq = y/2glo ~ 
1. 0733 X 10-2 m/s. 

Eq. ()26|) yields J\f relaxation times r = Tj and the same number of eigenvectors which 
evolve as 

N, = N,(0)expH/r,). (27) 
The number of particles in each discrete (gravitational) state Ni decreases as 

N, = J2N^3i^)^^Pi-t/^3)^ (28) 

3 

where Nij is the i-th component of the j-th eigenvector. In other words, we have to project 
the initial distribution of neutrons over the gravitational states, (A^^i (0) , N2 (0) , ...) on the 
eigenvectors Nj. Then the overall neutron count is 

N{h) = J2 N^3 (0) exp [-LjVoT, (h)] . (29) 

h3 

where the time is limited by the time of flight between the plates t = L^/Vq ~ 2 x 10^^ s. The 
dependence {h) acquires a distinct structure only if there is a pronounced hierarchy in the 
set of relaxation times Tj (h) and if there are pronounced changes in the longest relaxation 
times Tj (h) as a function of the separation between the plates. 

It is clear that we should take into account only the relaxation times which are comparable 
to or much bigger than the time of flight t = L^/Vq. For much shorter relaxation times we 
do not have to do any precise calculations and just disregard the components of the initial 
distribution which are decaying as the corresponding eigenvectors. 

As we will see below, these longest relaxation times Tj (h) , which correspond to the 
particles in the lowest energy states far away from the absorption edge, experience rather 
steep changes at certain values of the spacings between the plates hi < h2 < h^... However, 
this feature can be observed only for specific types of the roughness of the scatterer surface 
when the values hi, h2, h^--- are sufficiently far away from each other. One of the crucial 

14 



elements here is the relation between the times Tjj' and rj°-*. If Tjj' is shorter than rj^"*, then 
the loss of neutrons occur as a result of scattering-driven gradual increase in gravitational 
energy (slow diffusion of particles between the states j with a gradual increase in the state 
index j) until the neutron gets into the absorber. If this is the case, the disappearance times 
for neutrons from the neighboring states are comparable and it is not easy to extract the 
information on the properties of quantized gravitational states from the neutron count. 

In the opposite case of very short times Tj^\ the neutron gets absorbed essentially as 
a result of a single scattering act. In this case, one can disregard the interstate transitions 
and Tj ~ Tj^\ This regime is more advantageous for experiment especially if there is a 
distinct hierarchy of times Tj^\ The transition between these regimes is determined by 
the lateral size (correlation radius) of surface inhomogeneities and the asymptotic behavior 
of the correlation function. All this information is contained in the transition probabilities 
Wjj' {pj,Pj')- In what follows we will pay particular attention to the experimental conditions 
necessary for the observation of the proper regime. 

Note, that the transitions between the states are largely suppressed, Xj^j/ ^ Tjj, when 
the lateral size of the surface inhomogeneities on the scatterer R is much larger than the 
size of quantized states (in this case, R ^ Iq, or, in dimensionless units, r = i?//o ^ 1) 0. 
Since the momentum change 6/3 as a result of scattering by inhomogeneities of the size r is 
of the order of 6f3 ~ 1/r and is very small for large r, it seems that in this case the interstate 
transitions are allowed only when the energy gaps between the states are small. If this 
were correct, the population of the higher states would have remained constant no matter 
how long one would have waited (the gaps between the states increase with increasing state 
number j approximately proportionally to j). However, this is not completely true because 
even at r = R/Iq ^ 1 there are certain values of the spacings between the plates hj at which 
the transition channels j < — > J + 1 open spontaneously [7| leading to a drastic decrease in 
the relaxation times: the scattering driven momentum change 5/5 ~ 1/r is sufficient for the 
transition between the states j and j' when 

Pj - Pj+i ^ I (A,+i - A,) /V~e ^ 1/r (30) 
or, in normal variables, when the spacing H satisfies the equation 

1 [e,+i (H) - e, (H)] / ^/E^) ^ I, (H) /R. (31) 
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However, it is highly unhkely that this spontaneous opening of mode couphng channels 
affects the neutron count in the existing experimental setups. The observation of this effect 
requires drastically different experiments . To avoid dealing with this issue, we will assume 
that the size of inhomogeneities r is relatively small. 

In order to find the dynamics of the neutron count, we need information about the initial 
distribution of neutrons over the vertical velocities (over discrete energy states Nj {t = 0)). 
This distribution reflects the neutron distribution over the vertical velocities in front of 
the slit and is determined by the expansion of the classical wave functions in front of the 
slit over the gravitational states. This can, in principle, be done precisely, but only if one 
knows the exact distribution of particles over the vertical velocities in the initial beam. 
Unfortunately, at present the distribution of neutrons over the vertical velocities cannot be 
measured as accurately as the horizontal distribution. However, even if this distribution were 
available, the exact mathematical problem would still have been too complicated because of 
the complex structure of the experimental apparatus. 

The analysis of experiment, which has been performed in Ref. jj], indicates that the 
initial distribution is approximately uniform. 



where N (0) is the overall combined density of neutrons in all gravitational states and Z is the 
total number of such discrete states between the plates with < Uc- Note, that the uniform 
distribution ()32|) is a non-equilibrium one and should evolve in time as a result of scattering 
by surface inhomogeneities even in the absence of the neutron absorption. The equilibrium 
distribution, which does not evolve as a result of non-absorbing scattering irrespective of 
the form of collision operator, is given by equation 



Luckily for us, the exit neutron count depends mostly only on the initial population of the 
lowest states. The lifetime of the neutrons from the higher states is much shorter and they 
get absorbed well before getting to the counter. The velocities Pi for the lowest states can be 
considered as constants, jSi ^ jSo, and, if there is no singularity in the velocity distribution in 
the beam at zero angle, the distribution for the lowest states is close to uniform, Eq. 



iV, (0) = iV(0) /z, 



(32) 



Nj (0) = N (0) pj/ Zj^Pi ' (0) /Pj = const 



(33) 



i<Z 
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anyway. In this case, the difference between the distributions (jHHj) and ((221) is neghgible as 
well. 

The last remaining issue before solving the transport equations (j25|) is to relate the 
transition probabilities Wjj> to the surface roughness (to the correlation function of surface 
inhomogeneities) . 



III. SCATTERING PROBABILITIES: WEAK ROUGHNESS 



Roughness of (one of) the walls leads to scattering of neutrons, which, in turn, results in 
the transitions between the states, broadening of the lines, etc. In some intermediate energy 
range, when the lifetimes are already short but the spacings between levels are not yet large, 
this broadening can even effectively transform the discrete energy levels into a continuum. 
However, when the matrix elements are still relatively small, one can always separate the 
effects of scatteringon the energy spectrum from the transition probabilities that affect the 
particle transport 6| and, in our case, the neutron count. 



1 established methods 
and our earlier 
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. Since the geometry of the problem is 
,[3], the calculation should be modified 



We will start from the case of slight roughness for which there are we^ 
of calculation of the surface- driven scattering probabilities 0] (see also 
application to neutrons in a gravitational trap [l 
somewhat different from our previous applications 
though the general method remains the same. 

The rough scatterer corresponds to the barrier of the height Uc with position z = H + 
^ {x, y) , where ^ {x, y) is the random function which describes the surface roughness, instead 
of its "ideal" position at z = H in Figure 1. If we neglect small gravitational potential mgH 
in comparison with the barrier Uc (the difference in scales is 10^!), the roughness-driven 
"perturbation" becomes 

V (z, X, y) = Uce{z-H -i)- U^e (z-H)^ (x, y) 5 (z - H) f/,. (34) 

The transition probabilities are given by the squares of the matrix elements of this per- 
turbation averaged over the surface inhomogeneities, 

(35) 

The matrix elements should be calculated using the wave functions 

^! = i){x, y) ijj {z) 



w^jy (q,q') = (IV'jpgjwOf • 
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where ifj {x, y) are the properly normahzed plane waves and the functions ipj {z) are given 
by Eqs. (0) - (P- 

The matrix elements of the perturbation (j34p are (cf. Ref. 



Vjji = J exp {ix [p — p') + iy [q — q')) ^ {x) ipj {z) ipf {z) [mgH — Uc] 6 {z — H) dxds (36) 
= -e {p-p',q-q')U,i^,iH)ij,, (H). 

The scattering probabilities (jH^j) are then 

W,AP.T.p',q') = {\V,„,y,,\')^ = ap-p\q-q')U',^]iH)^}iH) (37) 

where C{p ~ p' i1 " l') ' the so-called power spectrum of inhomogeneities - is the Fourier 
image of the correlation function of surface inhomogeneities C, (x, y) , 



C {x, y) = {i{xi, yi)i{xi +x,yi + y)) = A J ^(xi, yi)^{xi + x,yi + y)dxidyi, (38) 

and A is the averaging area. Finally, assuming the homogeneity of the surface in ?/-direction, 
i.e., the lack of the dependence of ( {x, y) on y, we get 

(p, g; p\ q') = S{q- q') W,,, (p, p') , (39) 
(p, p') = C{P- P') {H) (H) . 

For computations, we will use the most common Gaussian correlation function, 

C (x) = f exp (-xV2i?') , C (P) = v^^'i?exp {-p'R^/2) , (40) 

where the amplitude £ and the correlation radius R characterize the averaged height and 
lateral size of surface inhomogeneities. Note, that though in many cases the real rough 
surface can be different from the Gaussian one there are reasons to believe that the 

n 

roughness in experiment [1| is close to Gaussian. The results for the relaxation time can also 
be sensitive to the type of the asymptotic decay of the power spectrum of inhomogeneities at 



large p ~\. Below we will choose such a regime for which the interpretation of experimental 
data would be the least sensitive to this uncertainty in parameters of the surface correlator. 
In our dimensionless variables, the transition probability ()39|) , ()40|) acquires the form 

w^,, = V2^ellWr exp [- (/?, - (3,,f /2\ ul^] (h) i,], (h) , (41) 
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where t] = i/lo and r = R/Iq. The condition of weak roughness is 

r^ = £//o<l,r. (42) 

In experiment, the potential barrier Uc is very high and can be considered infinite in 
calculations of the roughness- driven scattering probabilities for transitions between the lower 
gravitational states. In this case, Eq. (jHUI) should be replaced by 

Vn' = (P ~p',q- q') i^- (H) ^P'y (H) (43) 

and Eq. dH} - by 



2 „2 

/ ^ l I 7 " *-• X I 1 I I / J ■ / J ■ / I 

4m^ 
1 

^0 

where we introduce the following scale tq for transition times: 

1 



{Pj - Pf) = jr-^lWr exp [- (/3, - /3,0' rV2] {H) ^'f, (H) (44) 



The form 1)411) with Uc should be used only when calculating the transitions that involve the 
upper energy levels and the continuous spectrum above the threshold. 

For the gravitational states ((Zj), the transitional probabilities ()44|1 become 

V^n' fe- - Pr) = y^^n' = exp [- - /?,,)' rV2] (46) 

X a|af, [Ai' {h - Xj) - SBi' {h - X^)]^ [Ai' {h - Xy) - SBi' {h - Xf)^ 

with Aj given by Eq.®. 

For the square well states (fTHj) , the scattering probabilities acquire the " standard" form 



Wn' - Pf) = ^^^n' = ^~ ~ ''/^^ ^^^^ 

We also need the probability of transitions between the gravitational state j and the square 
well state j': 

Wjj' {Pj - Pj') = ^Wjj, (48) 
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All of the transition probabilities above describe transitions between the discrete states. 
In order to complete our calculation of the matrix of the relaxation times , we have 
to get the expressions for the transitions from the bound states into the continuous states 
1/tj'^\ To avoid extended calculations, we will first assume that instead of the continuous 
spectrum we are actually dealing with particles in a large box L ^ H with a set of very 
close levels fc„ = 7m/ L. Then for the transition probability from discrete states j in our well 
into the "continuous" spectrum inside the box L one should use Eq. (jH}, 



-(0) 



u 



9 exp 



An<e 



[1 + C0S2 {h^^ I (1 - U,l\n)\ 



(49) 



An = Mc + {-nnlQ/Lf . 



Switching back from summation over kn = irn/L to the integration over the continuous 
energy variable A, we get 



1 



(0) 



( v/^^ - Ve-«c-A)%V2 



(50) 



T^ToPj Jo VA — Wc 1 + (1 + Mc/A) COS^ (/ia/A + Uc) 

In principle, the above equations for the transition probabilities are sufficient for solving 
the transport equations from the previous subsection and for calculating the absorption rates 
for neutrons in the channel. 



IV. TOWARDS STRONGER ROUGHNESS 



It seems much more natural to design an experiment with scatterers with strong roughness 
rather than with the slight one: scattering by strong roughness ensures, almost automati- 
cally, a significant turn of the neutron velocity which results in direct absorption of a neutron 
by a plate. Therefore, neutrons form all states, except for the lowest gravitational ones, for 
which the wave function does not reach the scatterer, can be absorbed by the plates as a 
result of, essentially, a single scattering act. Then the count of the exiting neutrons will give 
directly, almost without any ambiguity, the number of particles in the lowest states that 
are not scattered. In this case, the dependence of the neutron count on the width of the 
gap between the plates should be the sharpest. If one is interested solely in identification 
of quantum states, one should definitely try to ensure that the surface roughness is strong. 
However, if one is interested in extracting the quantitative information from the experiment 
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about the energetics of the system, one should seriously consider using slight roughness for 
which the particle dynamics can be described very accurately as it is shown above. 

In other words, if one is interested solely in observing the quantization of neutrons, one 
should use a scatterer with strong roughness; if the purpose is the use of quantized states 
for ultralow energy calibration, one should use only rough systems which are subjects to a 
quantitative theory. Therefore, it is important to discuss the restrictions imposed by the 
conditions of slight roughness and the possibility of extending the limits of applicability 
of our theory in this particular case to stronger roughness. There are two types of strong 
roughness: high-amplitude and high-aperture roughness. As we will see, our theory can 
be extended to the high-aperture low-amplitude roughness, but not to the high-amplitude 
roughness. In the last Section we will argue that the high- amplitude roughness should not 
be used in experiment either. 

The surface roughness is considered to be strong if one of the conditions of slight rough- 
ness. 



breaks down. In such cases, one cannot usually say much except that one surface collision 
is sufficient for a complete dephasing of a particle. In our problem such an outcome could 
even have been looked at as beneficial since this would immediately allow us to estimate the 
corresponding relaxation time approximately as 



However, such an assertion is not always correct since the scattering probability W depends 
not only on the parameters of the surface (the correlation function of surface roughness () 
but also on the probability of the particles to get close to the wall which enters W via the 
particle wave functions and/or their derivatives on the wall (see, e.g., Eqs. (HT|) , (Hlj) for W 
for weak roughness). When such a probability is low, one can (and should) sometimes deal 
with the strong roughness as if it is weak. 

The calculation of the scattering probability consists of two steps: calculation of the 
matrix elements of the roughness-driven "perturbation" Vjj' and the calculation of the scat- 
tering probabilities W using these matrix elements. Since the latter step is easier to analyze, 
we will start from it. 



T] = i/lo < l,r 



(51) 



Tj = H/V,j. 



(52) 
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This second step is also less restrictive. Here one has to express the scattering probability 
W via the matrix elements Vjj>. The scattering probability is always given by the expression 
of the type (j35p , but with the scattering T-matrix T instead of the potential distortion V. 
The link between these two matrices is given by the operator equation 

f = V + TGV (53) 

where G is the Green's function. For quantized gravitational-like states j, j' (the lower 
states with h — Xjj' ^ 1), the position of the wall is deep under the barrier where the 



- (2/3) {h - Xf^ 



Then the 



wave functions are attenuating exponentially, roughly as exp 
matrix elements of V are exponentially small and Eq. ()53|) leads to T ~ resulting in the 
validation of Eq. (jH^j) even when the distortion V, by itself, is not small. This means that 
Eq. can be used as long as we are interested solely in the transitions to and from the 
gravitational states for which the matrix elements are always small even for not very small 
roughness. One should be much more careful with the transitions between the higher states. 

The siraplest way of evaluating the matrix elements Vjjr seems to be a generalization of 
approach 11 1 (see also the first Ref. ^) for the roughness-driven perturbation 

V (z, X, y) = U,e{z-H-0- U,e (z-H). (54) 
Instead of using Eq. ()34|) . we now write V as an expansion 

V (z, X, y) = U,Y. ^—r^^'"''^ - H) C {x, y) (55) 
1 ^■ 

with the matrix elements 

V,y (q - q') = f/, 5^ ^ / i,, (z) (z - H) (z) dz (56) 

I ^' Jo 

X j C (s) exp [i (q - q') • s] ds. 

The first integral reduces to a set of higher derivatives of the wave functions on the wall, 

or, in dimensionless variables, to /g ""''"^^/'j'"'' {h)tp^^ ^ ""^ {h). There- 
fore, if the amplitude of the corrugation is i, the expansion (|56|) becomes an expansion in 
i/lo- If the wave functions on the walls are (exponentially) small, as it is for the lowest grav- 
itational states, the matrix elements Vjj> can remain small even for not very small values 
of the corrugation amplitude £. However, in this case the all the constants in the expres- 
sions for the scattering probabilities change: instead of the wave functions on the walls one 
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should write a proper combinations of (h) from Eq. ()56|) and, more importantly, the 
correlation function of the surface corrugation should be replaced by the Fourier images of 
the higher moments of the correlation. 

If, on the other hand, the amplitude of the corrugation is small, i Iq, Eq. leads 
to exactly the same expression for the matrix elements as Eqs. fl3(i|l and, therefore, to the 
scattering probabilities (|35 p . (1411) . (1441) even when the aperture of the inhomogeneities is large. 
This means that, as long as the matrices V are small and are close to the T-matrix, V ^ T , 
the high-aperture roughness could still be treated in the same way as weak roughness while 
the high-amplitude roughness cannot. 



The above analysis on the basis of a generalized approach ll[ is not completely accurate. 
A a more accurate analysis can be done using the mapping transformation method (see 
both references iGll), which involves the lateral derivatives of the surface roughness in a more 
natural way than the method of Ref. This method involves the first lateral derivatives 
^i^y, which have the order of magnitude £/R and not I/Iq, explicitly. As a result, it looks 
that the expressions for the matrix elements Vjj> for the high-aperture roughness are different 
from Eq. . However, the collision operator in the transport equation involves not all the 
matrix elements of Vjj/ (q, q') , but, since the scattering probabilities are always accompanied 
by the energy 5-functions 6 {E — E') , only those that satisfy the energy conservation law. 
This means that Eq. ()35|) can be rewritten for our purposes as 

(q, q') ^ (9 + {W^v^rvW?), S (e, + e, - |-) . (57) 



It is relatively easy to see l&\ that the terms with the derivatives of ^ contribute only to 
the matrix elements with different energies E and E' if the amplitude of roughness is small, 
i <^ Iq. Therefore, the expressions for the matrix elements sX E = E' keep their weak 
roughness form even for the strong roughness with high aperture, but with low amplitude. 

Summarizing, the scattering probabilities remain the same as long as the amplitude of 
roughness rj = I/Iq is small, while the aperture of roughness rj/r = i/R can become large, 

< £ < /o- (58) 

This means that Eqs. (j41|) . (j44p for the scattering probabilities are valid not only for the weak 
roughness (jl^ , (jHU , but also for a much stronger high-aperture roughness , but only to 
the extent that one is interested in the transitions to and from the gravitational states with 
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exponentially small wave functions on the rough wall. In this case one can often disregard 
the correlation exponents exp [— [jSj — /5j')^r^/2] in the expressions for W and make the 
rest of the calculations as if the roughness were weak. Physically the disappearance of these 
exponents means that the changes of momenta in scattering are now unrestricted. However, 
the moment the transitions between the higher states come into play, these equations loose 
their accuracy. 

On the other hand, the large amplitude of roughness always leads to a completely different 
picture than the weak roughness. Note, that if we are interested in identification of the 
gravitational states when the distance between the walls H is comparable to the level size 
Zq, H > /q, the condition i <^ Iq is a necessary condition for the existence of the gravitational 
states. Therefore, we have to assume that this condition holds and that the only important 
type of strong roughness is the one given by Eq. (j3H|) . In the case of high-amplitude 
roughness with i 3> Iq, Eq. ()35|) still holds for the gravitational states and the structure 
of Eqs. (PT|) ■ ()44|) remains the same, though, instead of the correlation function of surface 
inhomogeneities ( [ttj — rcji), the expressions for W contain the summation over the higher 
order correlators. However, as it was explained above, we are not interested in this situation. 
What is more, under these conditions the state energies experience such large roughness- 
driven changes that the classification of states on the basis of an ideal well looses all meaning. 

For practical purposes, the results of this section mean that the our theory is still valid 
for a relatively strong roughness ()58p. but only in application to the lowest gravitational 
states and not to the square wall states. Since in this paper we were interested mostly in 
the direct transitions from the gravitational states into continuum, these results still hold 
for 1 ^ rj ^ r and the results in Figure 11 below are still within the applicability of the 
theory and the step-wise depletion of neutrons can be observed. 

V. RELAXATION TIMES 

A. Direct absorption vs. protracted diffusion between the states 

There are two ways how the neutrons can be absorbed by the plates (transit into the 
continuous spectrum) as a result scattering-driven interstate transitions. The first way 
is a relatively protracted diffusion between the states which results, eventually, in neutron 
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getting into continuum. The second way is by direct transitions from the initial state into the 
continuum without involving any intermediate states. For the latter method to dominate, the 
probabilities of direct transitions into continuum should dominate over transitions between 
individual discrete states. Obviously, this direct process has an undeniable theoretical appeal 
since in this case we do not need to solve the set of 120h x 120/i transport equations (|25p 
and only calculate the transition times for direct transitions into continuum Tj'^K 

However, as we will see below, following this route has important advantages for exper- 
iment as well. If the direct transitions dominate, the results will be the least sensitive to 
the correlation characteristics of the surface roughness (both to the correlation parameters 
and to the shape of the correlation function) and to the details of the initial distribution 
of neutrons over the energy levels. Then the interpretation of experimental data should be 
more reliable and unambiguous than is the case of protracted interstate diffusion. 

We will start from delineating conditions under which the direct transitions into contin- 
uum from all discrete states, including the lowest, are possible and dominant. The next step 
is the calculation of the corresponding transition rates which will immediately lead us to the 
results for the neutron count. 

The most relevant time scale is the time of flight 

t = L^/Vq ~ 2 X IQ-^s. (59) 



Numerically, the scale of the relaxation times in ()26|) is given, according to Eqs. 
by the coefficient tq ()45|) in the scattering probabilities VF, 

1 V2^ 



where t] = i/lo is the dimensionless amplitude of the surface corrugation, and t/ro ~ 2?)rf / [3. 
Since Iq ~ 5.871 /im and the spacing between the walls in experiment does not exceed 50 
/im, the reasonable limits on the amplitude of inhomogeneities are 10^'^ ^ V = ^/^o ^ 0.1, 
and 

2.3 X 10-^ < t/ro < 0.23. (61) 

(in existing experiments, 0.03 ^ ?7 ^ 0.1 and the lower limit in the inequality (jfjlj) is about 
0.02). 

The full relaxation times differ from Tq"^ by the dimensionless factors w. To simplify 
the problem, we have to follow the largest of these w. Note, that according to the above 
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expressions for W (see, e.g., Eq. (@7j)) the most probable transitions are the transitions into 
the highest allowed (by the correlation exponent) states. Therefore, the direct transitions 
into the continuous spectrum, when allowed, have the shortest transition times. What is 
more, the corresponding inverse times Tj'^^^^ contain summation over all the accessible states 
in the continuum and could become noticeably shorter than the transition times between 
the individual discrete states. 

There are two types of factors entering the dimensionless probabilities w: the correlation 
exponents with the negative indices — {(3j — /?j')^r^/2, which can only decrease w, and the 
factors related to the values of the wave functions on the walls. The correlation exponents, 
in general, encourage the transitions with the smallest change in momentum, i.e., the tran- 
sitions between the nearby states, unless, of course, the correlation length r is small (see 
below) and it does not matter. The coefficients originating from the wave functions on the 
wall, on the other hand, favor the direct transitions into (and from) the highest states (fac- 
tors j'"^ /h^ in Eqs. fl47|) , PH|) ). The balance of these opposing factors determines whether 
the attenuation of the gravitational states occurs via the gradual diffusion between the states 
or via direct transitions into the highest states. 

The possibility of the transitions is determined by the momentum transfer in the exponent 
of the correlation function, 

(/3,-/3,0'rV2= (v/^^- v/^^)%V2. (62) 

This means that the direct transitions from the lowest states into the continuous spectrum 
(and all transition between the discrete states) are not suppressed if 

- v^lTx^)' n,rV2 < 1, (63) 

(x = Uc/e) or, using the numeric value for Mc, 

r < 3. 8 X 10-=^Fo {x), R< 2- 2270 x lO'^Fo (x) y^^m, Fq (x) = v^/ (l - v^T^) (64) 

where the function Fq (x) is given in Figure 2. 

This condition is not overly restrictive. Unless the overall velocity of the beam is close 
to the absorption threshold, i.e., unless x = Uc/e is close to 1, the condition (j64j) can easily 
be satisfied. Essentially, condition (jMj) determines whether the particles, which are initially 
in the lower discrete states, can disappear as a result of direct transitions into continuous 
spectrum (or into short-lived upper states). 
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FIG. 2: Function Fq (x), Eq. 



Since the transitions between the lower states are much less likely than the transitions 
between the lower and higher states and between the higher states, we can use expressions 
PH|l for all Tjj} for the gravitational states and Eq. for all square well states. 

B. Transitions from the square well states into continuous spectrum (direct ab- 
sorption) 

We should start from estimating the times rj°^ for transitions into continuous spectrum 
fl50|) : if these times are short, at least for some discrete states j, t/rj'^^ ^ 1, this will help in 
evaluating the relaxation times for all other states as well. First, since Uc is large, one can 
replace cos^ (^h^/X + Uc) in the denominator of Eq. (j50|l by 1/2, 



(e — Xj) r^/2 > (e — Uc) r^/2 ~ lO^r^. If r is small, r < 10 (-R < 5 nm), as suggested in 




A)VV2 



(65) 



Since Xj < Uc < e, the index in the exponent changes from ( ^/e — A j — — u^) r^/2 to 
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FIG. 3: Function Fi (x), Eq. 



the previous section, the index in the exponent is always small and the correlation exponent 
can be disregarded, 

1 _u%rij]{h) r-^^ dX 2 _uj\ri;]{h)^ r^'^~^ z 1 

rf^ TrTo(3j Jq ^/XTu'c 3 + Mc/ a 7rro/3j Jq 3z + 1 y/1 + z 

(66) 

' ' ' r./N..^.. 1.2/v^- 1.76 + 0.245 In ^-1 



10 -^Fi ix) , F, (x) ^ 0.11 



3 - 2.45v^ 



The function Fi (x), is plotted in Figure 3. 

For the estimates, the function Fi (x) is mostly irrelevant except for e ^ Uc ^ 1^ 
According to Eq. (fTTjl . for deep square- well levels 

1 4.6 X lO^V 2A~ 



-F, (x) . (67) 



(for higher levels, one should use for ip"^ {h) a more accurate Eq. (HSJ). 

Taking into account Eqs. , (jUUI) for the time of flight t and characteristic time tq, it 
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is convenient to represent the ratio of the time of flight to the absorption time as 

' - ioVrSk(io%/«=)^iW- (68) 



rf ' h (10-3/3, 



Since even for the lowest square-well states Xj/uc > 10^^, this means that neutrons from 
all the square-well states, even for the lowest ones, disappear directly into the continuous 
spectrum during the time of flight at any reasonable amplitude of roughness and small 
lateral size of inhomogeneities as long as the correlation radius of surface roughness r is 
sufficiently small and the absorption threshold is not close to the overall kinetic energy (the 
ratio X = ""c/e is not close to 1). 

At larger size of inhomogeneities, the restrictions on the momentum transfer become more 
important and the integral in Eq. should be evaluated more carefully. For the lower 
states Xj <^ e, the absorption time is given by Eq. ()66p in which the integral Fi {x) should 
be replaced by the function F (x, ^/u^r) : 

1 uj^lor'^'^ (h) 

10 r ' F (x, v^r) , (69) 



r. 



(0) 



3 

rl/x-1 ^ ^ 

F(x,y) = 0.2 / dz exp 

^ ^ Jo 3z + iy^TT^ ^ 



The plots of the function F (x, i/u^r) at two values of r, r = 0.01; 0.001 are given in the 
Figure 4. The same function F (x, ^/^uZr) is plotted as a function of r in Figure 5 at 
X = 0.15; 0.3; 0.5. As it can be seen in the plots, choosing proper values for the correlation 
radius of inhomogeneities r and the ratio of the threshold energy to the overall kinetic energy 
is crucially important for having effective absorption. 

Function F (x, ^/u^r^ decreases rapidly with increasing x and even more rapidly with 
increasing r. Since y/u^ ~ 370 and the integrand depends exponentially on y/u^r, the 
dependence of F (x, y/u^r) on r is very steep: when r > 0.03 the momentum transfer re- 
strictions are so severe that the direct transitions from the lower levels to the continuous 
spectrum are suppressed by orders of magnitude unless x is very small [in existing experi- 
ment, X ^ 0.15]. For example, when r = 0.01 the function F (x, y/u^r) reaches 10"^ only at 
X ^ 0.26. 

With the same numbers as in Eq. ()68|) . the ratio of the time of flight to the time of direct 
transitions into the continuous spectrum is 

^ - ioV7:^^Sk(io%K)^(^'V^^)- (70) 
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FIG. 4: Function F (x, V^r), Eq. (jMJ, at three values of x, X = 0.15; 0.3; 0.5. 



C. Direct transitions from the gravitational states into the continuous spectrum 

Probabilities of direct transitions from the gravitational states into continuous spectrum 
are given by the same Eq. , but with different values of ipj {h), Eq. (fT^ instead of Eq. 
fll5|) . Since the integral remains the same, the only substantial change is the replacement of 
^ {h) in the r.h.s. of Eq. ^ by ^ , 

;M-10V|^&.i^(x,v/^r), (71) 

bj = W%B]/2 = lO^Oj [Ai' {h - Xj) - SjBi' {h - Xj)Y /u^ 
^ 0.3aj [Ai' {h - Xj) - SjBi' {h - X^)]^ 

where aj is given by Eq. (jSJ, 

a, = { (Ai' (-A,) - SBi' {-X,)Y - (Ai' {h - A,) - :SBi' {h - A,))^ . 

This means that the neutron in the gravitational states disappear on the time of flight 

2 

unless the coefficient aj [Ai' {h — Xj) — 5'jBi' {h — Aj)] and/or the function F (x, y^r) are 
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FIG. 5: Function F (x, ^r), Eq. JMl), at three values of x, X = 0.15; 0.3; 0.5. 



very small. This is so when the distance between the walls h significantly exceeds and 
the wave function on the wall and, therefore, bj are exponentially small. To evaluate this 
effect we need to calculate numerically the eigenvalues \j (h), Eq. ((Tj) and the coefficients 
bj, Eq. JH}. 

The dependence of the nine lowest (gravitational) energy levels on the spacing is presented 
in Figure 6. The crosses on the curves mark the points where the "size" of the level becomes 
equal to the distance between the walls (in our notations, X = h) and the wave functions 
start "touching" the walls or, in other words, the levels move from the gravitational to the 
square well domain. 

The coefficients bj (h) , Eq. (f7T|l . are presented in Figures 7,8 in two different scales. As 
we can see, the curves are extremely steep. This means that when scanning the distance 
between the walls h one encounters the critical values hj at which t/r^^^ changes from very 
small to very large values. According to numerical data, bj (h) changes between 10~^ and 
10-2 when 5.87 > h > 4.13 for bi, 6.64 > h > 5.76 for 63, 8.62 > h > 7.12 for 63, 
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FIG. 6: Nine lowest energy levels as a function of the spacing between the walls h. The stars 
mark the positions where the size of the energy of the level becomes equal to mgH (in dimensionless 
variables, = h) and the level shifts to the "square well" domain. 

10 > h > 8.34 for 64, etc. The change of b for all these levels from 10""^ to 10~^ requires 
change in h by only 0.2. Since such change in h is equivalent to the change in the distance 
between the walls by about 1.1 fim and the neutron count depends exponentially on the 
coefficients b, the neutron count should, under proper circumstances, become in experiment 
a truly step-wise function of the clearance between the walls. However, since the functions 
b {h) for higher levels overhang over each other, it is rather unlikely if within this observation 
technique one can resolve more than three or four gravitational levels. 

D. Neutron count 

As long as the main route for the neutron absorption is the direct transition into contin- 
uous spectrum, the neutron count at the exit is 




(72) 
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FIG. 7: The coefhcients bj (h) , Eq. (|71|) . which determine the square of the wave function on the 
surface of the scatterer, Bj. 



The count, of course, depends on the initial number of neutrons in each state, Nj (0). When 
the only long-lived states are the few lowest ones with the energies much lower than the 
overall kinetic energy of particles in the beam and than the scale of the energy dispersion, it 
is reasonable to assume that the initial populations of these states are the same, Nj (0) = No, 
and 



a function of the time of flight t, but as a function of the clearance between the walls h. 

Below we plot the functions n (h) and nj (h) for the overall neutron count and individual 
level depletions. According to the definition (f7^ . the neutron count from the individual 
level rij {h) always changes from 1 at large /i to at small h. The numerical value of the 




(73) 



while for each level 
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FIG. 8: The same as in Figure 7 but on a smaller scale. The stars on the curves mark the same 
points as in Figure 6. 

overall neutron count n (h) also provides the information about a number of non-depleted 
levels: for n {h) > 1, more than 1 level is still occupied; for n (h) > 2, more than 2 levels are 
still occupied; for n{h) > 2 - more than 3 levels are still occupied; etc. 

Since the functions hj (h) grow by orders of magnitude when h changes from 6 to 1, the 
ratios t/rj^^ also change by orders of magnitude and the neutron count n {h) should decrease 
very steeply. In principle, if the absorption times rj'^^ for different levels are vastly different, 
the curve n {h) (fTHjl should become a step-wise function with each step representing the 
decay of an individual state rij (6), Eq. (f7i|) . Unfortunately, this is not what happens. 

First, we will plot the function n {h) which describes the steepest achievable depletion of 
the neutron population, while still remaining within the applicability of the weak roughness 
theory. In other words, we will plot n (h) which corresponds to the largest combination 
of coefficients in Eq. ()70|) . Since the weak roughness theory is valid for roughness with 
the amplitude t] smaller than the correlation radius of inhomogeneities r, the largest t/rj^^ 
with given r is achieved at r] = r. The function F (x, y^r) , on the other hand increases 
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FIG. 9: Function lO'^r^F (x, for Eq. dZTJ) at three values of x, X = 0.15; 0.3; 0.5. 



with decreasing x though in experiment this parameter does not go below 0.15. The plot 
of the function lOV^F (x, y^r) , which enters the ratio t/Tj^\ Eq. (f70|) . as a coefficient, is 
presented in Figure 9 for three values of x- Clearly, the optimal conditions for experiment 
are x = 0.15, r] = r = 0.015 when lO^r^^F (x? v^^) ~ 0-23. 

The neutron count ()73|) with these values of x, ^ is presented in Figure 10 as a function 
of the interwall clearance. The same figure also gives the individual depletions of the nine 
lowest lowest quantum levels (jTlj). 

This figure clearly shows that though the depletion is sufficiently steep for an unambiguous 
demonstration of the quantization of the gravitational states (the total count is steep and 
goes to zero at /i 7^ 0). However, the "overhang" of the coefficients bj {h) [Figures 7,8] is 
too large for identification of the individual levels: though the individual depletion rates are 
clearly distinguishable, the overall neutron count is still a featureless function of h. This 
"overhang" can be compensated only by an increase in other coefficients in Eq. (f7n|) . This 
result means that the use of a scatterer with slight roughness is insufficient, even under the 
most favorable circumstances, for an observation of a step-wise change in neutron count 



35 



FIG. 10: The neutron count N/N (0), Eq. (|7.3)) . and individual depletions of the nine lowest energy 
levels, Eq. 1)741) . as a function of the spacing between the walls h at x = 0.15, rj = r = 0.01. 

and, therefore, for an identification of higher states, unless parameter x becomes noticeably 
smaller. 

On the other hand, the use of large roughness such as, for example, an increase in the 
amplitude of wall roughness to the value 1] = 30r, would modify the neutron count in Figure 
10 to the one depicted in Figure 11 with very clear steps that separate the depletions of the 
individual levels. Though the parameters themselves r = 0.015, rj = 0.45 are reasonable 
from the point of view of experiment, the large ratio rj/r = 30 is out of range of the slight 
roughness theory. However, the theory can be extended, as it is explained in Section IV, to 
a low-amplitude high-aperture roughness and can cover the situation plotted in Figure 11. 

Though the existing experiment does not produce as unambiguous picture as in Figure 
11, the similarities indicate that in experiments the roughness was not small. Additional 
information on the evolution of neutron count as a function oi rj/r can be found in Figure 
14 below. 
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FIG. 11: The neutron count N/N{0), Eq. H73|) . and individual depletions of the lowest energy 
levels, Eq. ()74() . as a function of the spacing between the walls h at x = 0.15, r = 0.015, rj = 30r. 
All higher levels are already depleted. 

E. Alternative geometries 

Similar experiments have been performed in alternative geometries with, for example, a 
mirror on the top and the rough scatterer on the bottom. The theoretical description of the 
change in neutron population with time remains the same as above with the only difference 
that all transition probabilities are now proportional to [^/^(O)!^ instead of \ip{H)f. For 
upper states this does not make any difference since the wave functions are symmetric and 
|'?/'(0)|^ = {ip {H)f . This is not so for the gravitational states (fT^ for which the wave 
functions on the top and on the bottom plates can differ by orders of magnitude. 

The ratio of the time of flight t to the relaxation times rj°^ for the gravitational states in 
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FIG. 12: Pairs of coefficients bj (h) (lower curves) and Cj {h) (upper curves), Eqs. (|7T|) . ((75|l . wliich 
describe tlie square of the wave function on the surfaces s = h and s = for the lowest nine states. 
Merging of these coefficients at small spacing between the walls h correspond to the transformation 
of the gravitational states into the symmetric square well states. 



this alternative geometry is given, instead of Eq. (fT^ , by the following expression: 



4^2 



r 



(0) 



10-3/3^ 
Ai' (-A) - 5Bi' (-A) 



(75) 



n 2 



Ai' {h-\)-S,B-^ (h-X,) 

As expected, the coefficients Cj are larger than the coefficients bj, especially when the spacing 
between the walls is large. The difference between these coefficients disappear when the 
spacing decreases and the states become more symmetric as it is clearly seen in Figures 12, 
13. 

When the coefficients Cj and bj are comparable, the corresponding depletion rates for 
both geometries are also of the same order of magnitude. Under the same conditions as in 



38 



FIG. 13: The same coefficients bj [h) and Cj [h) , Eqs. (|7T|) . ((75|) . for tlie lowest four levels as in 
Figure 12, but on a smaller scale, which describe the square of the wave function on the surfaces 
s = h and s = 0. Merging of these coefficients at small spacing between the walls h indicate the 
transformation of the gravitational states into the symmetric square well states. 

Figure 9, though the depletion rates for the first level in direct and inverse geometries are 
different because of a noticeable difference in ci and bi, the depletion rates that include nine 
levels do not differ much between themselves (Figure 14; a pair of curves marked d — 1 and 
^-1). 

However, with an increase in coefficient ()75p. the difference in depletion rates between 
"direct" and "inverse" geometries, which is due to the differences between cj and bj, becomes 
more and more pronounced. Figure 14 contains several more pairs of curves for direct 
and inverse geometries which correspond to the amplitude of roughness equal to r]/r = 
2; 3; 4; 5; 6; 8 (the curves are marked accordingly). For even larger roughness, such as in 
Figure 11, the depletion rates for the inverse geometry is so fast (the scale is 10~^^!) that 
the curves for direct and inverse geometries cannot be plotted in the same figure; even for 
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FIG. 14: The depletion rates for direct {d; the scatterer is on the top) and inverse (i; the scatterer 
is on the bottom) geometries for x = 0.15, r = 0.015. The curves are marked by the ratio rj/r = 
1; 2; 3; 4; 5; 6; 8 and by the geometry, d or i. 

rj = 8r the numbers for inverse geometry are already too small. Note, that when the spacing 
between the walls becomes small and the energy levels go sufficiently high up into the square 
well domain, the pairs of curves, which describe the depletion rates in direct and inverse 
geometries, merge between themselves. 

Figure 14 illustrates the evolution of neutron count with increasing aperture of roughness 
rj/r. According to this figure, the increase of the ratio rj/r from 1 to 8 is still insufficient to 
produce a step-wise plot of the neutron depletion with well pronounced steps as in Figure 
11. However, at ?7 = 8r there are already some signs of the appearance of the first step on 
the depletion curve. 
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VI. PERSPECTIVES AND CONCLUSIONS 

A. What to expect from strong roughness experiments? 

On the face of it, the case for shifting experiments towards strong roughness is clear: the 
curves should become sharper and the chances of identifying the individual gravitational 
states higher. However, as usual, there are the pitfalls. 

The strong roughness can come in two flavors: large-amplitude roughness, comparable to 
the spacing between the walls or to the size of the gravitational levels, and the roughness 
with small lateral size of inhomogeneities in comparison with their amplitude (the high- 
aperture roughness). The use of the former one is rather dangerous. The main problem is 
that the distance between the walls becomes an ill-defined variable and, because all energy 
parameters are very sensitive to this spacing, the energies extracted from such experiments 
are unreliable. Therefore, though the sharpness of the curves improves, the curves, by 
themselves, can be used for nothing more than qualitative identification of the individual 

n 

states. The only exception might be the so-called adiabatic roughness [15|] in which the 
profile of the walls changes so slowly that the wave functions can adjust to the local shape of 
the walls. Still, even in this case the uncertainty in energies is large making the quantitative 
results useless. 

The low-amplitude high-aperture roughness, on the other hand, can both improve the 
resolution for the lowest states and allow one to extract useful data on the interaction 
between neutrons and solid surfaces. This can be achieved only when the dominant depletion 
regime is the direct scattering-driven absorption of neutrons. The moment the intermediate 
square-well states become involved, the accuracy of data extracted from experiments with 
high-aperture walls is lost. 

B. Angular resolution 

Another way of increasing the sensitivity of experiment to individual levels is to decrease 
parameter x = Uc/E which describes the ratio of the absorption threshold to the overall 
kinetic energy of particles: obviously, the smaller is this ratio, the higher are the chances 
of a particle to get absorbed after scattering. Mathematically this effect is described by a 
steep function F (x, ^/u^r) , Eq. , in Figure 4. 
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The experimental control over the value of x could be limited. However, there might be 
an easier alternative to lowering x- experiment, and, as a result, in the theory above, 
very little attention has been paid to the third spatial direction, y. What we were looking 
at was the scattering-driven change of direction of particles from the beam direction x to 
the vertical direction z which could provide the particle with sufficient kinetic energy mv1/2 
to overcome the absorption threshold Uc- Obviously, as a result of scattering the particle 
acquires not only the vertical velocity Vz-, but also the component of velocity Vy in the third 
spatial direction. Though this component of velocity does not affect the chances of the 
particle to be absorbed by one of the horizontal walls, the particle can, nevertheless, escape 
the exit counter if the velocity Vy and, therefore, the deviation from the initial trajectory, 
are large enough. On the other hand, the exit neutron counter cannot distinguish whether 
the neutrons have escaped as a result of absorption or because of angular dispersion of the 
beam in ^/-direction. Therefore, instead of decreasing parameter x, one can achieve the same 
effect simply by decreasing the size of the counter in ^/-direction. Of course, the price of 
using the low-angle counter will be the overall drop in the number of detected neutrons. 
Therefore, following this route might require an increase in beam intensity. 

From a theoretical standpoint, the use of the angular resolution requires taking into ac- 
count the scattering in y direction instead of averaging in Eq. ()21|). Since we have the full 
2D expressions for the scattering probabilities W , this will lead just to computational com- 
plications rather than to conceptual difficulties. For a wide aperture counter this extension 
would hardly change the nature of the results. However, for narrower counters, as it is ex- 
plained above, this might turn out to be important. Essentially, we have to introduce a new 
angular counting threshold and find the proper replacement for the function F (x, ^/u^r) . 
The structure of the calculations will remain the same especially if we limit ourselves, at 
least in the beginning, to the calculation of direct high angle scattering processes rather 
than looking at protracted diffusion of particles over the angles. 

C. Contribution from the intermediate states 

In this paper we investigated only the situation when the direct scattering into continuum 
dominates the neutron depletion. In principle, one should take into account the processes 
that include the indirect absorption processes that involve the intermediate energy states. 
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When the direct absorption processes are efficient, namely, when r is sufficiently small, 
adding diffusion over intermediate states should not speed up the absorption. Involving 
the intermediate states can be beneficial only for some combinations of parameters when 
the direct processes are suppressed because of a relatively large value of r. The absorption 
via the intermediate states can sometimes compensate for the presence of the correlation 
exponents with Ucv"^ that slow down the direct absorption processes. On one hand, the 
probability of scattering-driven transition Aj Xj is proportional to the square of the 
wave function on the wall which decreases rapidly with lowering of the state energy Xj 
and, therefore, favors transitions with as large energy increase as possible. On the other 



nearby states. The competition of these opposing trends determines at what values of r the 
intermediate states are important for absorption. 

The depletion rates for higher "square well" states are always very high. Therefore, at a 
relatively large value of r there could even be a range of parameters in which the depletion of 
the gravitational levels via the intermediate state could become much more efficient than the 
direct transitions over the threshold that the value of the threshold could become irrelevant 
and the neutron count could lose its dependence on x- 

In this paper we were focused on the search for "optimal" roughness that ensures the 
steepest neutron depletion when the wave functions start "touching" the walls; under these 
conditions, the multistage processes via intermediate states do not seem to be necessary. 
However, the quantitative description of experiment cannot be complete without a more 
careful look at the role of intermediate states. We plan to address this issue later in a 
separate paper. 

D. Short-range forces near the surface 

A fundamental issue is whether it is possible to extract the parameters of the (ultra)short- 
range forces acting on neutrons near the solid surfaces (and whether such forces exist at all). 
Our approach to particle diffusion along rough walls, in which the surface roughness is treated 
in the same way as bulk perturbations, is well-suited to answer the question. The short-range 
forces near the surface can be considered as another 5— type perturbation in addition to the 
roughness- driven perturbation ()H4|1 . The matrix elements of this additional perturbation 



hand, the correlation exponent 




favors transitions into 
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are determined by the same equations as the ones that we used for the calculation of the 
surface-driven transition probabilities. The regular part of these matrix elements determines 
the shift of the eigenstates A and the values of the wave functions on the walls, 6j, Cj which 
are necessary for the calculation of the neutron depletion rates. After such modification 
the calculation can, in principle, proceed along exactly the same lines as it is done above. 
Whether or not this program will be implemented in the future will be determined by 
an experimental success in observation of a well-pronounced step-wise dependence of the 
neutron count on the spacing between the walls. If this high resolution is not achieved, the 
inclusion of the corrections from the short-range forces would be meaningless. 

E. Conclusions 

In summary, we studied the possibility of detection of the quantized gravitational states 
of neutrons. We applied our theory of quantum transport along rough surfaces to a neutron 
beam between flat and rough walls and calculated the roughness-driven transition probabil- 
ities between the quantum states. This allowed us to found the depletion rates due to direct 
absorption process as a function of the spacing between the walls and parameters of the wall 
roughness. From the theoretical standpoint, the main achievement is the extension of our 
theory to systems with strong, high-aperture roughness and the application of the theory to 
systems with absorbing, "leaky" walls. 

Prom the point of view of experimental applications, the focus was on flnding the optimal 
roughness parameters that ensure reliable resolution of several gravitational states. The main 
stumbling block is the "overhang" of the wave functions from different states which, in most 
situations, prevents resolution of the nearby states. Apart from the roughness parameters, 
the neutron count is also very sensitive to the ratio of the particle kinetic energy to the 
absorption threshold. The most suitable experimental conditions are those in which the 
characteristic times rj^^ for the direct transitions into continuous spectrum (absorption) are 
short in comparison to the time of flight except for the lowest gravitational states for which 
the probability to flnd the particle near the scatterer /absorber are exponentially small. 

According to our results, the preferable conditions for the observation of the well pro- 
nounced gravitational states are the following: 

• Weak roughness is sufficient for estabhshing the fact of quantization, but is not suffi- 
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cient for resolving the individual quantum states. The sharpest experimental results 
within the weak roughness regime are expected when the amplitude and the correlation 
radius of weak roughness are £ ~ ~ 0.08 //m. 

• The resolution of quantum states requires the use of stronger roughness. The use of 
low-amplitude high-aperture roughness Iq, H ^ i > R is preferable to high amplitude 
roughness i > Iq or i > H. In the latter case, though the levels can be resolved, it is 
virtually impossible to extract useful quantitative information from the experimental 
data. 

• According to our results, resolving two or three lowest states would require i > 20R 
when using low-amplitude high-aperture roughness. 

• In the case of high aperture roughness with rj ^ r one can even neglect the correlation 
exponents and make the rest of the calculations as if the roughness is weak. Physically 
the disappearance of these exponents means that the changes of momenta in scattering 
are now unrestricted and the results are now not very sensitive to the correlation radius. 

• The threshold velocity for the absorption should be noticeably lower than the overall 
beam velocity, x = Uc/E < 1. This condition, if necessary could be replaced by the 
condition that the width of the detector in y-direction Ly (or the width of the plates) 
is considerably smaller than the distance to the detector (length of the plates) L. In 
this case, the absorption energy Uc is replaced by the threshold for the disappearance 
of the neutrons in y-direction, Ey ~ {Ly/Lf E . 

We also formulated several suggestions for future theoretical and experimental work which 
could help in reliable identification of the gravitational states and use the acquired informa- 
tion for the study of the short-range forces near solid surfaces. 
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APPENDIX A: FIGURE CAPTIONS 

Figure 1. (Color online) Schematic drawing of the potential well in vertical direction 
z. Dashed lines - the "real" potential; solid lines - the potential used in the paper. Since 
mgH/Uc ~ 10^^ and we are particularly interested in the lowest energy levels, the corrections 
are negligible [with an accurately scaled potential, the difference between the dashed and 
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solid levels cannot be seen]. 
Figure 2. Function Fq {x), Eq. 
Figure 3. Function Fi (x), Eq. (pj)) . 

Figure 4. (Color online) Function y^r), Eq. (j69j) . at two values of r, r = 

0.01; 0.001. 

Figure 5. (Color online) Function y^r), Eq. at three values of x? X = 

0.15; 0.3; 0.5. 

Figure 6. (Color online) Nine lowest energy levels Aj as a function of the spacing between 
the walls h. The stars mark the positions where the size of the energy of the level becomes 
equal to mgH (in dimensionless variables, Aj = h) and the level shifts to the "square well" 
domain. 

Figure 7. (Color online) The coefficients bj (h) , Eq. (f7T|) . which determine the square of 
the wave function on the surface of the scatterer, 

Figure 8. (Color online) The same as in Figure 7 but on a smaller scale. The stars on 
the curves mark the same points as in Figure 6. 

Figure 9. (Color online) Function lO^r^F (x, ^/vZr) for Eq. (|7T|l at three values of %, 
X = 0.15; 0.3; 0.5. 

Figure 10. (Color online) The neutron count N/N (0), Eq. ((7Hj). and individual depletions 
of the nine lowest energy levels, Eq. ()74p . as a function of the spacing between the walls h 
at X = 0.15, r] = r = 0.01. 

Figure 11. (Color online) The neutron count N/N (0), Eq. (f7H|) . and individual depletions 
of the lowest energy levels, Eq. (ffljl . as a function of the spacing between the walls h at 
X = 0.15, r = 0.015, 7] = 30r. All higher levels are already depleted. 

Figure 12. (Color online) Pairs of coefficients bj (h) (lower curves) and Cj {h) (upper 
curves), Eqs. (f7T|). (fTSj) . which describe the square of the wave function on the surfaces 
s = h and s = for the lowest nine states. Merging of these coefficients at small spacing 
between the walls h correspond to the transformation of the gravitational states into the 
symmetric square well states. 

Figure 13. (Color online) The same coefficients bj {h) and Cj {h) , Eqs. (f7T|) . (f73j) . for 
the lowest four levels as in Figure 12, but on a smaller scale, which describe the square of 
the wave function on the surfaces s = h and s = 0. Merging of these coefficients at small 
spacing between the walls h indicate the transformation of the gravitational states into the 
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symmetric square well states. 

Figure 14. (Color online) The depletion rates for direct {d; the scatterer is on the top) 
and inverse (i; the scatterer is on the bottom) geometries for x = 0-15, r = 0.015. The 

curves are marked by the ratio rj/r = 1; 2; 3; 4; 5; 6; 8 and by the geometry, d or i. 
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